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Abstract 

We prove that a topologically generic network (an open and dense set of networks) of three or more 
inhibitory neurons have periodic behavior with a finite number of limit cycles that persist under small 
perturbations of the structure of the network. 

The network is modeled by the Poincare transformation which is piecewise continuous and locally 
contractive on a compact region B of a finite dimensional manifold, with the separation property: it 
transforms homeomorphically the different continuity pieces of B into pairwise disjoint sets. 

PACS 2008 codes: 87.19 lj, 87.19 11, 87.19 lm 
MSC 2000 92B20, 34C25, 37G15, 24C28 
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1 Introduction 

We study the dynamics of an abstract dynamical system modeling a network composed with n > 3 neurons, 
where n is arbitrarily large, that are reciprocally coupled by inhibitory synapsis. 

Each neuron is modeled as a pacemaker of the type integrate and fire [?] , including the case of relaxation 
oscillators. The internal variable Vi — Vi(t) describing each neuron's potential, for i — l,2...,n evolves 
increasingly on time t during the interspike intervals: Vi(t) has positive first derivative dVi(t)/dt > and 
negative second derivative d 2 Vi(t)/dt 2 < being the solution of a deterministic autonomous differential 
equation of a wide general type. 

When the potential Vi reaches a given threshold value, the neuron i produces an spike, its potential Vi 
is reseted to zero. It is suppose that all the neurons are inhibitory, there are no delays and the network is 
totally connected. When the neuron i spikes, not only its potential Vi changes, being reseted to zero, but 
also, through the synaptical connections, an action potential makes the other n — 1 neurons j / i, suddenly 
change their respective potentials Vj with a jump of amplitude — H{j < 0. 

The instants of spiking are defined by the evolution of the system itself, and not predetermined by 
regular intervals of observation of the system. We analyze the state of the system immediately after each 
spike, in the sequence of instants of spiking of the network. We prove that the state of the system after 
each spike is a function F of the state after the prior spike. This function F is the so called Poincare map. 
The study of the dynamics by iteration of the Poincare map is not an artificial discretization of the real 
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time dynamics. On the contrary, the dynamics of F and its properties (for instance periodicity, chaotic 
attractors) are equivalent to those of the system evolving in real time. 

In |MS-1990] the technique of the first return Poincare map to a section transversal to the flux was first 
applied to study neuron networks, in that case, an homogeneous network of excitatory coupled pacemakers 
neurons. In [BCRG-1996] the same technique is applied to networks of inhibitory cells, analyzing the real 
time dynamics via a discrete Poincare map F. This map is locally contractive and picccwise continuous in 
a compact set of R n_1 . We include the proof of these properties in the section 2 of this paper. 

In Section [4l we prove Theorem 14. li which is the main abstract mathematical result: 

Locally contractive piecewise continuous maps with the separation property generically have only periodic 
asymptotic behavior, with up to a finite number of limit cycles that are persistent under small perturbations 
of the map. Generic systems have a topological meaning in this paper: they include an open and dense 
family of systems. 

As a consequence we obtain the following applied result: 

Generic neuron networks composed by n > 3 inhibitory cells exhibit only periodic behavior with a finite 
number of limit cycles that are persistent under small perturbations of the set of parameter values. 

This is a result generalizing the conclusions obtained for two neurons networks in |BTCE-199i] and 
|UB-1992j . 

On the other hand non generic dynamics are structurable unstable: they are destroyed if the system is 
perturbed, even if the perturbation is arbitrarily small. We refer to those as bifurcating systems. 

The results of this paper are proved in an abstract and theoretical context, using the classical qualitative 
mathematical tools of the Topological Dynamical Systems Theory. The systems have discontinuities and 
evolve in finite but large dimension n. It is still mostly unknown, the dynamics of discontinuous and large 
dimensional systems. That is why in sections 3 and 4 we include an abstract theory of topological dynamical 
systems with discontinuities. The piecewise continuity and the local contractiveness help us to obtain the 
thesis of persistent periodicity in the Lemma l4.2l of this paper. On the other hand, the separation property 
will play a fundamental role to obtain the thesis of density, and genericity of the periodic behavior in the 
Theorem 14.11 

Even being our mathematical analysis theoretically abstract, we observe that the conclusions about the 
dynamics of our model of neuron networks, fit with those obtained by experiments in computer simulations 
with mutually coupled identical neurons in networks of up to 10 10 cells, as reported in the following papers: 

In (PVB-2008] it was observed the transition among different periodic activity, indicating that the sim- 
ulation data perfectly fit to experimental and clinical observations. In the computer simulated experiments 
the alterations of the discharge patterns when passing from one periodic cycle to another, arise from changes 
of the network parameters, changes in the connectivity between cells, and also of external modulation. 

In |PWVB-2007j the computer simulated experiment shows the dynamics of the network of a large 
number of coupled neurons. It was observed to be significantly different from the original dynamics of the 
individual cells: the system can be driven through different synchronization states. 

Our thesis of Theorem 14.11 is only applicable to deterministic systems. Nevertheless their conclusions 
also qualitatively fit with computer simulations of neural systems with randomness [IV-2007], |TV-2007) . 
which also show the generation of detectable preferred firing sequences. 

2 A mathematical model of the inhibitory neurons network. 

We include the detailed proof of the mathematical translation from a physical model of n inhibitory 
pacemaker neurons network to the dynamics of iterations of a piecewise continuous contractive map F : 
B i— > B, locally contractive and with the separation property, as first posed in |BCRG-1996| . The model is 
applicable for any finite number n > 3 of neurons in the network. 

The phase space of the system is the compact cube Q = [—1,1]™ C R™. A point in the phase space is 
V = (Vi, V-j, . . . , V n ), describing the potential Vi of each of the neurons i 6 {1, 2, . . . , n}. We assume that 
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the the phase space is normalized: the threshold level of each of the neurons potentials is 1, so 1 is the 
maximum of Vi . Also the minimum Vi is normalized to — 1 , and that the reset value of Vi , after a spike of 
the neuron i, is 0. 

Definition 2.1 The physical model. The point V in the phase space Q evolves on time t, during 
the interspike intervals of time, according to an autonomous differential equation and changes without 
delay in a discontinuous fashion in the exact spiking instants, according to a reseting-synaptical rule. The 
two regimes, during the interspike interval, and in the spiking instants respectively, are precisely defined 
according to the following assumptions: 

12.11 1 Inter-spike regime assumptions. Vi(t) is the solution of a differential equation 

^ = 7<(V i ), 7*:[-M]~R> 7,eC\ n(Vi)>0, Yi(Vi)<0 Wi £[-!,!}. (1) 

where C 1 denotes the space of real functions in [—1,1], continuous and derivable with continuous deriva- 
tive in [—1,1]. 

The assumption 7; > reflects that each neuron potential in the inter-spike interval is strictly increasing 
while it does not receive interactions from the other neurons of the network. This comes from the hypothesis 
that each isolated neuron i is of pacemaker type, i.e. from any initial state 1^(0) S [—1, 1), the potential 
spontaneously reaches the threshold level 1 for some time t = tj > 0, if none inhibitory synapsis is received 
in the time interval [0, U]. 

The assumption 7^ < 0, which we call the dissipative hypothesis reflects that the cynetic energy E c = 
(l/2)(dVi/dt) 2 is decreasing on time while the potential freely evolves during the interspike intervals. In 
fact: dEjdt = (dV/dt)(d 2 V/dt 2 ) = 7z(^M^M^) < 0. 

The most used example of this type of inter-spike evolution is the relaxation oscillator model of a 
pacemaker neuron, for which ji(Vi) — —ctiVi + [3i where < cm < are constants. For this type of cells 
the differential equation {T]) is linear, and its solution can be explicitly written: 

V t {t) = (A/a,) - m/ai) - Vi(0)]exp(-ait). 

12.11 2 Consequences of the inter-spike regime assumptions. We define the flux 

&(y) = (<S>{{V 1 ),,¥ 2 (V 2 )...,<S> t n {V n )) 

as the solution with initial state V — (Vi, V2, ...,V n ) of the differential equations system given by |T]). 
Precisely: 

^ffill= 7i ($*(^)) Vt, ^(Vi) = Vi (2) 

As 7i £ C 1 we can apply the general theory of differential equations to deduce the following results, as a 
consequence of the assumptions in |T|) : 

• Two different orbits by the flux do not intersect. 

• If B and A are two (n — l)-dimensional topological and connected sub-manifolds of E™ transversal to 
the vector field 7 = (71, . . . , 7„), then the flux transforms homeomorphically any set of initial states 
in B onto its image set of final states in A. 

• For each constant time t it holds the Louville formula: 



= exp\ i^KV^ds. (3) 
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12.11 3 Spiking instants computations. For each initial state V G Q the first spiking instant t(V) in 
the network is defined as the first positive time such that at least one of the neurons of the network reaches 
the threshold level 1. This means that 

t(V) = min U(Vi), where = 1 t = (4) 

l<i<n 

J(V) = {ie{l,2,...,n}:t(V)=t i (V i )} (5) 

is the set of neurons that reach the threshold level simultaneously at the instant t(V). It is standard to 
prove that for an open and dense set of initial states there is a single neuron i reaching the threshold level 
first, i.e. #J(V) = 1, J(V) = {i}. 

12.11 4 Spiking-synaptical assumptions. In the spiking instant t the reseting and inhibitory synaptical 
interaction without delay produces an instantaneous discontinuity 

a : <P J ^{V) ^<j(<f t{v) {V)) 
in the state of the system, according to the following formulae: 

• If #J(V) = 1, {i} = J(V) then = 1 and: 

ai(<f>l (V) (Vi)) = (spiking-reseting rule) (6) 
a) (<$>} V) (Vj ) ) = max {-1, ®J V) (Vj) - H^} Vj ± i (synaptic rule) (7) 

where Hij > is constant, depending only on and gives the instantaneous discontinuity jump in 
the potential of neuron j ^ i produced through the inhibitory synaptical connection from neuron i 
to neuron j 

• If #J(V) = k > 2, {ii,i2, ■ ■ ■ ,ik} — J{V) then a is multiply defined, having k possible vectorial 
values a 11 , a 12 , . . . , a lk , where a lh is defined according to formulae ((6]) and ([7|). 

We also assume that the network whose nodes are the cells and whose sides are the synaptical inhibitory 
interactions Hij , is a complete bidirectionally connected graph. Precisely: 

< e = min (8) 

12.11 5 Relative large dissipative assumption. We assume the following relations between the 
functional parameters 7.; in the differential equations ([T]) governing the dissipative interspike regime, and 
the real parameters Hij in the formula ([7]) governing the spiking-synaptical regime. 

maxtfy < \ (9) 

I <*IA\ Ci I AW / min ' min V,£ [1/4,3/4] \iSVi)\ n „v 

max |7i(3/4)-7 i (3/4)| < l -j — ^ (10) 

max^j- H zj _ 1 < minj min y ie [1/4,3/4] |7~(^)| 
minijtj Hij 4max, 7* (3/4) 

Condition @ assumes that the discontinuity synaptical jumps Hij are not relatively as large as the 
widest range [0, 1] of the potential of the cells when they act as oscillators between the reset value and 
the threshold level 1, free of synpatical interactions. 
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The hypothesis (jTOJ) and pip verify for instance for homogeneous networks in which all the functions 
7i and all the synaptic interactions Hij are constant independent of the neurons But as they are open 
conditions, they also verify if the network is not homogeneous but the neurons and the synaptical jumps 
are not very different. Finally they also verify for networks that are very heterogeneous, but the dissipative 
parameter of the system mini m i n v i e[i/4 ,3/4] l7i(^) is large enough. 

The assumptions above can be also possed for some number < a < 1/2 instead of 1/4 in inequality 
([9]), the number 1 — a instead of 3/4 in the values of Vi, and 2/(1 — 2a) instead of the denominator 4, of 
inequalities (fT0|) and (fTTj) . Nevertheless, and without loss of generality, in the computations of this work 
we will take the assumptions above with a = 1/4 to fix the numerical bounds. 

2.2 Comments about the physical model. 

The hypothesis ([9]), (fTOf and (|11| will allow to prove the so called separation property in Theorem 12.91 
in this paper. This property is essential to prove that the family of all the systems which exhibit a limit 
set formed only by a finite number of limit cycles is dense, which leads to the topological genericity of such 
systems. 

We observe that the assumptions in {lj , (JB]) and Q are more general that what they a priori seem. In 
fact, if instead of the variables Vj which describe the electric potentials of each of the neurons, we used 
other equivalent variables, the vector field 7 of the differential equation fT]). and the synaptical vectorial 
interaction a given by (O and 0, would have other coordinate expressions. 

For instance, each isolated cell % acts as an oscilator, whose potential Vi varies in the interval [0,1]. 
We can diffeomorphically change the variable Vi to a new one V £ [0, 1], called the phase of the oscilator, 
which by definition, evolves linearly with the time t, during a time constant Tj. In the new variables the 
differential equation governing the phase state Vi will be dVi/dt = 1/tj and the flux will be linear in Q. 

In [BCRG-1996] it is developed the model in such phase variables Vi for which the flux is linear, and 
it is defined the synaptical inhibitory interaction jumps — Sij < in the phase state Vj, when the phase 
Vi reaches the threshold level 1. To be equivalent to the constant jumps —Hij in the old variables Vj, it 
is showed in BCRG-1996 that the interaction jumps — Sjj < in the new phase variables Vj, must be 
functions Si t j(Vj), strictly increasing with Vj and such that Vj — sy(Vj) is also strictly increasing. In a 
widest model the functions Sij(Vj) are continuous but not necessarily differentiable. 

In resume, up to a change of variables, the model assumed in this paper in hypothesis fl}, (J6]) and 
includes for instance the model in [BCRG-1996 in which the flux is linear during the interspike interval 
regime, and the synaptic jumps in the spiking instants adequately depend of the phase of the postsynaptic 
neuron. 

Definition 2.3 The Mathematical Model. In this subsection we will define a Poincare section B C Q 
of the dynamical system modeling physically the network of n inhibitory neurons defined in 12.11 We then 
shall define the first return Poincare map F : B t—> B. We will prove that this map is piecewise continuous, 
locally contractive and has the separation property. These properties justify the Definition 12. 15[ at the end 
of this section, in which we will model and analyze this kind of inhibitory neuron networks through the 
abstract mathematical discrete dynamical system defined by the iterates of its Poincare map F. 

12.31 1. The Poincare section B. Let B G Q = [—1, 1]™ be the compact (n — l)-dimensional set 
defined as follows: 

n 

B = (J B k where B k = {V e Q : V k = 0} (12) 

k=l 

The topology in B is defined in each B k as the induced by its inclusion in the (n — 1) dimensional subspace 
{Vk = 0} of R n . Each B k is transversal to the flux defined in 12.11 2 solution of the system of differen- 
tial equations (JTJ) , because the vector field 7 in the second term of this differential equations has all its 
components strictly positive. 
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After each spike, the state of the system is in B, due to the reset rule in equality ©. So the system 
returns infinitely many times to B from any initial state V *E Q. 

The geometric illustration of the three neurons system. In figure Q] the cube Q is represented 
for n = 3 dimensions. The three coordinate axis Vi, V2, V3 are the three edges of the cube that are hidden 
in dotted lines, at the rear part of the cube, intersecting pairwise orthogonally in the unique hidden vertix 
O of the cube. The axis of Vi goes from O to the front, the axis of V% to upwards, and the axis of V3 to 
the right. 

The Poincare section B is formed by the three faces of the cube that correspond to at least one of the 
potentials Vi, V2, V3 equal to zero. B is the union of the three faces of the cube at the rear part, that would 
not be seen if the cube were not transparent, each one in a plane, orthogonal in the origin, to the respective 
coordinate axis. 

The increasing orbits of the flux inside the cube, are in that figure, parallel lines orthogonal to the plane 
of the figure. Due to the perspective each of the orbits is seen as a black dot in the figure. Each black dot, 
for instance "a", represents a linear segment of an orbit. The dot "a" is the orbit of the flux $*(V) from 
the initial state V = (Vi, V>, 0) e B, in the left face of the cube (that would not be seen if the cube were 
not transparent), to the front side of the cube, where the neuron i = 1 reaches the threshold level 1. 

In that moment, the neuron i = 1, whose potential increased to reach one (the state of the system is 
in the front vertical face of the cube), resets to zero, and the state of the system goes (through a dotted 
horizontal line parallel to the axis Vi), from the vertical front face of the cube (where V\ = 1) to the parallel 
vertical rear face in the Poincare section B, where Vi = 0. (in the figure: from the point "a" at front, to 
the point "b" at back). 

In that vertical rear face belonging to B, the other two neurons j = 2, 3, whose potentials V2, V3 were not 
reset, suffer a reduction of their potentials, of amplitudes H12 and -H13 respectively, due to the inhibitory 
synaptic rule. That is why, the system does not stay in the point "b" of the figure, but jumps to "c" , 
always in the rear face of the cube, corresponding to V\ = 0. 

At that instant, immediately after the first spike, from the point "c" in the backward rear face of the 
cube, the system starts to evolve again according to the differential equation, in the inter-spike regime, 
moving on an orbit inside the cube. In the figure this orbit corresponds to a segment orthogonal to the 
plane of the observer, collapsed in the black dot "c" due to the perspective. This new orbit arrives to the 
upper face of the cube, (also in the black dot "c" of the figure), meaning that neuron i — 2 arrived to the 
threshold level one. 

One could believe that figure [T] is too particular, because the flux is linear inside the cube Q, with 
orbits that are parallel lines. But, due to the Tubular Flux Theorem, any flux tangent to a vector field 
7 = (71,72,73) such that 7, > 0, after an adequate differentiable (generally non linear) change of variables 
£ in the space, becomes a linear flux whose orbits are parallel lines, and project orthogonally to a certain 
plane. One could transform all the dynamical system in these new coordinates £(V), but in this case the 
synaptic jumps Hij in the old variables Vj would become synaptic jumps sy in the new variables £j(V). 
Maybe the matrix (sy) becomes dependent of the new variables V\ = £i(V), Vi — ^(V), ... ,V n = £„ (V), 
similarly to what was remarked in 12.21 

The Tubular Flux Theorem is also valid in any dimension n > 2, and this geometric model has all the 
data that we will develop analytically in this paper. In particular, the Tubular Flux Theorem and the 
linearization of the orbits, are used and analytically written in the proof of Theorem 12.111 

The observer of the systems does not need to "see" the orbits inside the 3-dimensional cube Q of the 
figured] to study its dynamics. He or she just see a point jumping inside the plane hexagon on which the 
cube Q projects orthogonally to the flux, and orthogonally to the plane of the draw, from the viewpoint 
of the observer. There is a transformation F from this hexagon to itself, giving the position of the point 
"c" from the initial state "a" , then the position of "e" from "c" , etc. The dynamics by iterates of this 
transformation F describes exactly the same dynamics of the system, which indeed evolves with continuous 
time t, but is disguised as discrete. It is not the system which was discrete nor the observer who made it 
discrete. The observer positioned in the adequate viewpoint, without modifying the system, just to see it 
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in its discrete disguise. 

Nevertheless there is a problem with this discrete system, if two or more neurons got the threshold level 
simultaneously. For instance if the point "a" would be in one of the three frontal sledges of the cube, that 
are sledges inside the hexagon, marked as not dotted lines in figure[U there would be more than one possible 
consequent state F(a) = "c", depending on which frontal face the system chooses to reset and apply the 
synaptic rule. That is why the transformation F has discontinuities, and is not a uniquely defined map in 
the frontal sledges, or lines of discontinuities. 

All the arguments in this section could be obtained geometrically in the n-dimensional cube Q, just 
after its projection on an adequate poligon on an hyperplane, in which the state of the system evolves 
accordingly to a discrete transformation F. 

Figure 1: Model of a 3 neurons network in the 3-dimensional cube: Reaching the threshold level of neurons 1,2 and 3 
corresponds to the front faces 1, 2 and 3 respectively of the cube. Points marked in black correspond to the linear evolution 
from backward faces to the front faces. Firing of neurons 1, 2 or 3 correspond to the jumping dotted lines from the respective 
front face to its parallel backwards face. The figure shows the evolution after 8 spikes of the neuronal system: 1,2,1,1,1,1,1 
and 3. 




12.31 2. The partition of B in the continuity pieces Bi. Recalling the definition of the spiking 
instant t(V) in equalities ([4]), and the definition of the set J(V) of all the neurons that reach the threshold 
level at time t(V), in equality ([5|), we define the following subset Bi of the Poincare section B, for any 
i = 1, 2 . . . , n: 

Bi = {V E B : iE J(V)} = {V E B : t(V) = U(Vi)} (13) 

In other words, the set Bi is formed by all the initial states V in the Poincare section B such that the 
neuron i reaches the threshold level before or at the same instant than all the other neurons of the network, 
from the initial state V. 

From the implicit equation at right of formulae we deduce that Bi is compact, and that its interior 
int(Bi) is formed by all the initial states for which ij(Vi) < tj(Vj) for all j ^ i. Then int(Bi) f]int(Bj) = 
Vi ^ j. 

As the flux is strictly increasing inside Q, from any initial state V E B there exists a finite time 
t(V) defined by equalities ([4]). Therefore V E Bi for some not necessarily unique i E {1,2,..., n}. Then 
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B = U"=i Bi an d the family of subsets {Bi}f =1 is a topological finite partition of B (i.e. it is a covering of 
B with a finite number of compact sets whose interiors are pairwise disjoint.) 
The compact sets Bi are called continuity pieces. 

We define the separation line S of the partition {-Bi}™ =1 , or line of discontinuitiesas the union of the 
topological frontiers dBi of its subsets Bi. Precisely: 

n / n \ 

S=\JdB l= |J (f?j nB j )=B\ i\J intBi (14) 

i— 1 ineqj \i— 1 / 

12.31 3. The first return Poincare map F. 

The first return map F : B ^ B to the Poincare section £? = [J™ =1 -Bi is the finite collection of maps 
fi : Bi^ B defined as 

fi(V) =o- l {& {v \V)) VVeB t . 

where <& is the solution flux defined in 12. H 2 of the system of differential equations (JTJ) , t(V) is the spiking 
instant defined by equalities (j4|) and a 1 is the synaptic vectorial map defined in 12.11 4. 

For simplicity we denote -F|b< = fi and, when it is previously clear that V £ £>;, we simply 
denote F to refer to the uniquely well defined map 

We observe that F is uniquely defined in 1J™ =1 intBi, and multi-defined in the separation line S. 
Applying the formulae Q, ([6]) and 0, we deduce: 

^kCV^/i^O = ((/<)i,(/i)a, •••.(/<)«) VFGA where 
(^|B i )iW = (/i)iW = - max{-l, 

(fk)i(V) = - max {-1, S* ^) - Hij} Vj (15) 

where by convenience we agree to define f/jj = +1. 

12.31 4. Piecewise continuity of the Poincare map F. 

The formula (|15p implies that /, = F|s 4 : Bi^ B is continuous, and, as Bi is compact, then fi(Bi) is 
also compact. 

The formula (TT5|) changes when one passes from £?i to £?/j with i ^ h, so F is multidefined in the points 
of S = {Ji-£ h (Bi n Bh). Besides F may be discontinuous in V° £ Bi n £?h because 



is not necessarily equal to 



lim F(V) = fi(V°) 



lim F(V) = f h (V°) = 

V£mtB h ,V^V a 



Remark 12.31 4: We agree to define the image set F(V) of a point V £ B as {fi(V) : i such that V £ 
Bi}. The image set F(V) is a single point if V £ int(Bi) because intBi does not intersect Bj for j ^ i. 
The image set F(A) of a set A C B is by definition = V) VeA F(V). 

12.31 5. The positive reduced Poincare section B + . We will define a subset B + C B such that 
F P (B) C £> + for all p large enough. Our aim is to study the limit set of the orbits, therefore the last 
property allows us to restrict F to B + . 

B+ = {V £ B : < Vi < 1 - e Mi = 1, 2, . . . , n} (16) 

where eo > is the minimum of the absolute values of the synaptic interactions Hij > for i ^ j, as 
assumed in 12.11 4. equality JHJ). Also, by hypothesis (J9j) we have 

0<e o <i, ^<i_ £o <i (17) 
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12.31 6. Properties of the positive reduced Poincare section B + . 

The set B + is homeomorphic to a compact ball in R"" 1 , (property that the whole Poincare section B 
does not have). In fact, B + is the union of n compact squares B~£ = {V G R" : Vk = 0, < Vi < 1 — eo} = 
{0} x [0, 1 - eo]" -1 C R" -1 such that, for h k: B+ n B% ^ is formed only by the (n - 2)-dimensional 
lines {Vk = Vh = 0} in the frontiers of both squares. 

The dynamics properties of F restricted to the positive Poincare section B + , which precisely justify the 
restriction to B + , will be stated and proved in Theorems 12. 4[ |2~TI 12.91 and 12. Ill They justify the definition 
of the abstract mathematical model in 12.151 whose dynamics in the future and attractors will be studied 
in the following sections. 

12.31 7. The positive continuity pieces B^ of the Poincare map. We define 



where B + is the positive reduced Poincare section defined in 12.31 5 and Bi are the continuity pieces of the 
Poincare map F, defined in 12.31 2. equality (fl3)) . and iri 12.31 4. 

Remark: It is rather technical to prove that Bf is homeomorphic to a compact ball in R n_1 . We 
sketch here a proof, leaving the technical details: Bf is the pre-image in B + by the flux <I> = $*(')(•), 
of the (n — 1) square Ai = {V £ Q : Vi = 1}. The flux is injective and continuous from B + onto its 
image $(i? + ) C A = [J A ., because it is transversal to A and to B + and two different orbits of the flux 
do not intersect. Any two orbits that intersect B + in different points do intersect A in different points. 
Continuous and injective maps <& from a (homeomorphic) ball B + in R"^ 1 onto a set $(f? + ) C R n , are 
homeomorphisms, due to the Theorem of the Invariance of the Domain. Therefore $(_B + ) is homeomorphic 
to a (n — l)-dimensional compact ball, contained in A. 

Also $(_B + ) (~l Ai is. To prove this last assertion, be aware that the intersection of two compact 
homeomorphic balls is not necessarily a single homeomorphic ball, but it holds in our model, because 
the frontiers of $(i? + ) and Ai have some symmetric properties due to the fact that the flux $ has n 
components each one depending only on the respective single variable Vk- 

Then, Bf is the homeomorphic image by <i> _1 of the compact (homeomorphic) ball $(i? + ) n Aj. □ 

Theorem 2.4 . The return map to the positive Poincare section B + . 

The positive reduced Poincare section B + C B defined in \2.3\ 5, is forward invariant by the Poincare 
map F\ B + : B + \— * B + , and it is reached from any initial state in B. Even more, 



such that, if V G B verifies Vi < 1 — eo V 1 < i < n, then the interspike interval t(V) > T . 

Proof: According to the formula f4|: t(V) = min^ U(Vi) where t = U(Vi) is the solution of the implicit 
equation Q\{Vi) — 1. We integrate the differential equation (pj with initial condition Vi, and recall that 
liiVi) > 0, while the real solution $f(Vi) < 1 is strictly increasing with s (for Vi constant) and it is 
the solution of an autonomous differential equation. Using the hypothesis V. < 1 — eo, and applying the 
inequality (fT7|) . we obtain: 



Bf = B+ n Bi 




Lemma 2.5 There exists a constant positive minimum time T 



maxfe 7fc(3/4) 




(18) 
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l = l-eo + / 7i ($f(l-e ))ds 

JO 

where < UiV) < U(Vi) and (V) = 1 - e 0) being U(Vi) the time that takes the flux $|(Vi) to be 

equal to 1 — eo from the initial state Vj < 1 — eo. 

Recall that 7i(Vj) > is strictly decreasing with V%. 

7i ($J(l - e )) < 7<(*i(l - eo)) = 7<(1 - eo) < 7,(3/4) Vs > 

rU(V t )-U(V t ) 

l<l-e + / 7i(3/4)ds 
Jo 

eo < 7,(3/4) MVi)-U{Vi)] < 7,(3/4) ^) 

U(Vi)>T= £ ° j— -t Vi, t(V)=maxti(V5)>r. □ 

maxfe 7fe(3/4) 

Proof of Theorem 12.41 It is enough to prove the following two assertions: 
Assertion [23 A: Fj(F) < 1 - e VV G B (even if V £ B+), V j = 1, 2 . . . , n. 

Assertion 12. 41 B: There exists a constant t\ > such that for all V G -Bi, if i^(V) < for some j ^ i, 
then i^(V) - V,- > ei. 



Note that the assertion 1 2. 41 B states its thesis in particular if V ^ _B , and also if V G -B + and V} = 0. 

Recall that from the formulae (fT5)) of the Poincare map F: Fi(V) = for all V E Bi. Observe that, 
being Vj > —1 for all V E B, from the Assertion [2]4jB we deduce that the first number p > 1 of iterates of 
F such that F P (B) C B + is at most equal to 1 + Integer-Part(l/ei). 



To prove the Assertion 12.41 A. apply the formulae (fl5|) of the return Poincare map F, and recall the 
assumptions ©, ©. If V G Bi then 

Fi(V0 = 0, F i (l/)=max{-l, sf^V}) - Hy } < 1 - min JJ« = 1 - e (19) 

To prove the Assertion 12. 4I B. fix V G £?; such that, for some j j^i 

Fj(V)<0 (20) 

Use the formulae (fl9|) . We assert that 

$< (y) (V}) < J (21) 

In fact, if it were greater or larger than 1/4, as ifjj < 1/4 due to hypothesis (J9j> , the formulae (119|) would 
imply that Fj(Vj) > contradicting our hypothesis (120| . 

Due to the hypothesis of the differential equation (fTJ), the function jj(Vj) is strictly decreasing with Vj, 
and the flux is strictly increasing with t. Use the integrate expression (TPSt of the differential equation, 

to compute $j\Vj), the inequality (f2"Tj) and the Lemma [2~5l to deduce: 

0< S <t(U) =* Vj<$ S j(Vj)<$f v \Vj)<\ =* 7i(^(^))>7i(V4) 
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,h ' V) 7i(«5(^))j da > 7i (l/4) > min 7fe (l/4) • T = 

Recalling the integral equation (fT8)) and the formula (jT9|) of the return map F, we deduce: 

Fj(V) - Vi > ^(Vj) - V S - Ha = f iV) 7iW(V^))j ds - ff« 



= ei 



, T/ s T/ / mint, 7fc (1/4) max,^ g tJ \ 

" Vj ~ £ ° Uax fc 7,(3/4) ^ J 

To end the proof it is enough to show that e\ > 0. Recall from equality (J8j) that eo = min^- Hj > 

ei _ min fe 7fe(l/4) max^iJ^ _ min fc 7fc(l/4) - max fe 7fc(3/4) / max 4 ^ iJ. t . 



eo maxfc 7^(3/4) minify ffy maxfc 7^(3/4) \ min^- H i: 

ei = 7 / t (l/4) - 7h(3/4) + 7h (3/4) - 7fc(3/4) _ /mnx,,., //,, _ 1 
eo max fe 7 fe (3/4) \ mill,, , //,, 

where we have taken h and fc such that 7^(1/4) = minfc 7fc(l/4), 7fc(3/4) = maxfc 7fc(3/4). 

ei ^ 7 h (l/4) - 7h (3/4) - | 7h (3/4) - 7fc(3/4)| /max,^ iJ y 



e maxfc7fc(3/4) \ mil1 '- i H ij 

ei > 7h(l/4) - 7h(3/4) - max^t |7fe(3/4) - 7fc(3/4) _ / max^ 3 Hg _ ^ 



e maxfc7fc(3/4) V. :nin,, , //, 

Applying the mean value theorem of the derivative of 7 /j, which is negative due to the dissipation hypothesis 
of the differential equation in assumption ([1]), we obtain 

3 1^\ , UM mini miny je[ i/4 )3 /4] |7i(^)l 



< 7ft (l/4) - ^(3/4) ={j-\y (-7'Jx))l xe[ i/4,3/ 



2 

The last inequalities and the assumptions (flQ|) and (|11|) of relative large dissipative, imply: 

fi minj min^g^/4,3/4] h^Q4)J f^_l_^\ _ Q D 
e maxfc 7 fc(3/4) ^2 4 A J 



Remark 2.6 Formula of the Poincare map in B + . 

As a consequence of Theorem 12. 4[ from now on we will restrict the Poincare map F to the positive 
section B + . In fact, from the statements of Theorem 12.41 it is deduced that the forward dynamics and the 
limit set of the orbits to the future, of the restricted F, will be the same as those of F in the whole Poincare 
section. 

Due to TheoremHUif V e B+ then F(V) C B+ . Therefore {F\ Bi )j{V) > Vi, j. Using the formula 
(|15p we can rewrite the expression of the Poincare map, without the maximum: 

F\ Bt (V) =fi(V) = ((fi)i,(fih,...,(fi) n ) WGB+ where 
(F\ Bt )i(V) = (fMV) = = ^ (Vi \Vi) - H u 

{F\ Bt ) J {V) = {h) J {V) = ^(V^-Hj Vj (22) 

Theorem 2.7 Local injectiveness of the Poincare map. 

The Poincare map F defined in formulae h2S\) , restricted to each of its positive continuity pieces Bf 
defined in 12.31 7. is injective. 
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Proof: Fix a continuity piece Bf of F in the positive Poincare section B + . The piece Bf will remain 
fixed along this proof. Therefore we will denote F instead of /j. 

Take V, W G B+.such that F(V) = .F(W) £ £+. We must prove that V = W. 

Due to the formulas (f22|) of the Poincare map: Fi(y) = Fi(W) = and 

= $f v \ Vj ) + H ti = zf^W) + ff« = V j 

As ify is constant, we deduce that 

Therefore the vectorial flux <&'(F) <E Q = [—1,1]™ defines an orbit from the initial V G B + that 
intersects the orbit from the initial state W G -B + . Two different orbits of the flux do not intersect. Then, 
the two orbits are the same. If necessary changing the roles of V and W, we deduce that 

$*o (V) = W for some < > 

B + C B, so V has at least one component Vk — and all of them not negative. But ^ is the strictly 
increasing in time solution of the differential equation d^/dt — 7^ (<&'■), with jj > for all j. 

We deduce that if V G B+, and if <f> to (V) = W for t > 0, then Wj > V j, and therefore W ^ B+. 
As we know that W £ B+ and t > 0, we conclude that t = 0, and then W = (V) = = V. □ 

Definition 2.8 The Separation Property. We say that F verifies the separation property if 

/,-(£+) n /,(!?+) = Vi^j 

where {S^}, i — 1, 2, . . . , n, are the continuity pieces of i 7, in the positive Poincare section B + , as defined 
in !2.3l 7.. and fi is the continuous expression of F\ B + according to the formulae (|2"2"]) . 

Note that Bf is compact for all i, and F is continuous in each Bf . Therefore the image F{Bf) is a 
compact set. Then, the separation property implies that there exists a minimum positive distance a > 
between the images by F of two different continuities pieces. 

Theorem 2.9 The Poincare map F verifies the separation property. 

Proof: Take Bf and Bf with i 7^ j. The formulae (|2"2"]) of the Poincare map F\ B + : B + cBn5 and 
the Theorem 12.41 imply that 

VVGB+: (/ i ) i (F)=0, (/i)j(V)>0 Vj^i 
V W G B+ : (/, ), (V) = 0, (/, MVO >0V^j 
Then /,(£?+) fU (£+) = 0. □ 

Remark 2.10 Global injectiveness of the Poincare map. 

From Theorems 12.71 and 12.91 it is deduced that the Poincare map F is globally injective in B + . In 
fact, if V 7^ W are in the same continuity piece Bf , then fi(V) ^ fi(W) because F\ B + — fi is injective. 

And iiV^W respectively belong to two different continuity pieces Bf and Bf for i ^ j, then fi(V) 7^ 
fj(W) because fi(Bf) n fj(Bf) = 0, due to the separation property. We deduce that if V ^ W then 
F(V) P| F(W) = 0, where the image set F(V) of a point is defined in the Remark [2~3l 4. 

Theorem 2.11 Local contractiveness. 

The Poincare map F\ B + is uniformly contractive, but not infinitely contractive, in each of its continuity 

pieces Bf . Precisely, there exist two constant real numbers < a < A < 1 and a distance dist in the 
positive Poincare section B + = [J Bf , such that, for all i — 1,2, ... , n: 

a dist {V,W)< dist {fi(V),fi(W)) < A dist (V, W) VV,W E Bf 

where fi : Bf B + is the continuous restriction of F to Bf , according with formulae S22\) . 
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Remark: The distance dist of Theorem 12.111 induces the same topology in B + as homemorphic to 
a compact ball of R™ -1 C R™. In fact, along the proof of the Theorem 12.111 we will construct a linear 
projection ir : R™ h-> R n_1 and a diffeomorphism £ : K™ i— > R™ of C 1 class, such that: 

dist (V, V + dV) = \\n(d£dV)\\ where || • || is a norm in R". 



Proof of the Theorem 12.111 The continuity piece B^is fixed. For simplicity of the notation, along this 
proof we will use simply F to denote fi. 

The existence of the distance dist and the contraction rate A is proved in the Theorem 3 of C-2008]. 
For a seek of completeness we include here some pieces of the proof of |C-2008j , adding to them the existence 
of the bound contraction rate 0<<7<1,<7<A. 

Due to the Tubular Flux Theorem there exists a C 1 diffeomorphism which is a spatial change of variables 
£ : V i— > V from Q C R n onto Q C R n , such that £,\b+ — id and the solutions of the differential equation 
© in Q verify 

dV/dt = a 

in Q, where a G R™ is a constant vector with positive components. It verifies: 

£W))=£(V) + <r-t, d£- 7 (V) = a VVgQ 

Define in R ra the ortogonal projection ir onto the (n — l)-dimensional subspace 

a x Vi + a 2 V 2 + ■■■ + a n V n = 

The flux of the differential equation (fT]), after the change £ of variables in the space, is ortogonal to that 
subspace, and is transversal to £(£> + ) = B + = B + (recall that is the identity map). 
Consider any real function g : R^™ -1 ) R: 

VV,V + dV€R n : w(dV) = 7r(dV + g(V)-a). 
VV,V + dV, U G B k , define dist (V, V + dV) — ||tt(^ dV)\\ 

dist {V, U) = [ ||7r(d£ v +t(ir-v) ■ (U - V)\\ dt (23) 
Jo 

It is left to prove that fi : Bf i— » B + is contractive with this distance. 
Let us apply fi to V and V + dV in e Bf . We use the equalities (|22|) . 

We shall use the Liouville derivation formula of the flux of the differential equation respect to its initial 
state: 

d&j/dVi = exp ( f * 7 j dsj) 

Define: 

— a = max max 7,-(Vi) < 0, — a* = min min 7-,-(Vi) <0 

3 Vjel-1,1] J 3 VS€t-l,l] 3 

Use the Lemma \2. 5 1 to bound uniformly above zero the inter-spike intervals t(V): 

< T < t(V) 

Recall that t(V) is the solution of the C 1 implicit equation $*^ V *J(T^) = 1. Then t(V) is a continuous real 
function of V E B + , and B + is a compact set. So, ?(V) is also upper bounded by a constant: 

t(V) < T* 
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Derive the formulae (|2"2")1 to obtain: 

F(V) - F(V + dV) =dF-dV= [(dFj/dVAdVj + {dF 3 / dV^dV^^ 

dFj/dV, = (d^W/dV^^ =e X p(/* W 7^(^(VS))d S )) 6 [e~ a ' T ' , e~° T } 

(24) 

fli^/aVi = (d^(V,)/dt)\ t=J(v) ■ {dU{V i )/dV l )=g{V)- l3 ^} v \v j )) 

where g(V) — dti(Vi)/dVi is the real function obtained deriving respect to Vi the implicit equation given 
in 1 = $**^*^(Vi). Call ej to the j— th. vector of the canonic base in 1" and join all the results above: 

7T • d£ (F(V + dV) - F(V)) = 7T ■ d£dF ■ dV = 
= K-d£ (EUidFj/dVi) ■ dVjej) + tt • <%(g(V) ■ 7 ($^>(V)) = 

= 7T • d£ (EU^Fj/dVi) ■ dV j e ) + g{V)ir ■ d£ ■ 7 (*^(V))) = 

= * • ^ (ff^dFj/dVj) ■ dV j e 3 ) + 7r(<7(F) • S) = tt • (E?=i(^/^) ' ^ ^) ( 2 5) 

We define the numbers a and A: < a = eT a T < e~ aT = A < 1 and observe from the computations in 
([24f that: 

< a = e- a " T * < dFj/dVj < e" QT = A < 1 
Applying the definition of the differential distance dist in (j2"3"| . and the equality (f2l>l) , we obtain: 

dist {F(V),F(V + dV)) = \\Tr(d£-dF-dV)\\<\\\Trd£-dV\\=X dist (V, V + dV) = A \\ir(d£ dV) || 
dist F(V + c2F)) = ||7r(d£ • • dV)|| > a ||7rd£ • dV|| = a dist (V, V + dV) = a \\ir{d£, dV)\\ 

Integrating by formula (j!?3"|) we conclude: 

a dist (V,Z7)< dist (F(V), F(U)) < A dist (V,C7) □ 

Remark 2.12 Local homeomorphic property of the Poincare map. 

Each continuity piece fi of the Poincare map in Bf is an homeomorphism onto its image. 

It is an immediate consequence of Theorem 12.111 and the global injectiveness of F. Even more, the 
continuous restriction fi = F\ B +, is Lipschitz with constant A < 1 and its inverse (defined from fi(Bi) i— > Bi) 
is also Lipschitz with constant l/er > 1. Then fi is an homeomorphism onto its image. □ 

We note that the same result can be obtained without using the Theorem 12.111 Due to Theorem 12. 71 
fi is injective, and due to the formulae (|22p . fi is continuous. Due to the Theorem of the invariance of the 
Domain, any continuous and injective function from a ball in R n_1 to K ra_1 is an homeomorphism onto its 
image. 

In the following corollary we resume all the conclusions of this section: 

Corollary 2.13 If the network of n inhibitory neurons verifies the assumptions of the physical model, 
evolving with real time t in the phase space Q C R™ as stated in (JTJ) , © , 10 , © , © , (fTU)) and (TIT]) , then 
there exists a Poincare section B + that is homeomorphic to a (n — l)-dimensional compact ball, and a 
return map F : B + i— > B + , which has the following properties: 

a) F is piecewise continuous. Precisely: there exist a finite partition {i? 2 + }i<i<n of the Poincare section 
B + , formed by compact sets B^ homeormorphic to compact balls of K n_1 , with pairwise disjoint interiors, 
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and there exist n continuous maps ft : ^> B + being F(V) = {/<(V) : i such that V £ Bf } for all 
V e £? + . 

As a consequence F is univoquely defined as fi in the interior of its continuity piece Bi, and multidefined 
as ft, ft in Bi n 5,, if i ^ j, / 4 ^ /_,-. 

b) F is locally uniformly contractive and not infinitely contractive, i.e. for some metric dist in B + 
the exist constants < a < A < 1 such that for all 1 < i < n: a dist (V,W) < dist (fi(V), ft(W)) < 
A dist (V,W) VV,W E B+. 

c) F has the separation property, i.e. ft(B^) n ft(Bj) = if i ^ j. Therefore, there exists < a = 
min^ dist (ft(B+),f j (B+)). 

Note that from b) and c), it is deduced that F is globally injective in B + , as proved in Remark 12.101 
Also from b) it is deduced that ft : B+ i-> ft(Bf) C B+ is an homeomorphism onto its image. 

2.14 Comments about the mathematical model. 

Due to Corollary |2.13i all the general results that we will prove for abstract piecewise continuous maps F 
verifying a), b), c), are applicable to the networks of inhibitory neurons in the assumptions of the physical 
model stated in 12.11 Nevertheless the reciprocal of the Corollary 12.131 does not hold. Given a map F 
verifying a), b), c) there does not necessarily exist a network of inhibitory neurons in the hypothesis of the 
physical model stated in 12.11 for which F is its first return Poincare map. 

Nevertheless we can wide our scenario of possible inhibitory neuronal networks models. In fact, the 
properties a), b) c) are open (in the uniform C° + Lipschitz topology of the finite family of maps ft). 
Thus they are not only verified by systems for which the differential equations ([1]) are independent in the n 
variables Vi, but also if the system is of the form dV/dt = j*(V), where 7* : K n 1— > R™ is a C 1 vector field, 
near enough the given 7 = (7 = 71, 72, . . . , j n ), even if 7* does not verify all the hypothesis stated in |T|). 

Also the matrix {Hi j)i j of synaptic interactions in the network can be substituted for any matrix 
(H*j)ij(V), not necessarily constant, but functions near the constant matrix (flijOij an d so, still verifying 
the assumptions ©, (fTU|) . (fTTj) . Therefore, without changing the synaptical rules in equations © and 
([7]), but allowing the synaptic interactions slightly depend of the postsynaptic potentials, we will obtain a 
Poincare map F still verifying the thesis a), b), c) of the Corollary 1 2. 131 

Besides, as observed in the subsection 12. 2[ the physical model includes looser hypothesis than those 
specified in 12. H modulus any differentiable change of the variables of the system. So, also in those models 
the properties a), b), c) are verified by an open family of systems. 

Finally, the properties a) b) c) of the Corollary 12 . 1 31 are verified by many other models, in which the 
interspike regime is stated as a dynamical system depending continuously on time t and on the initial state 
V , but not necessarily as regular as to verify a differential equation. The dynamics of the potential Vi in 
the inter-spike interval may be given by a flux defined continuously in time t, strictly increasing on 

t, continuous but not necessarily differentiable respect to t nor to the initial state. But not all such general 
models are in the aim of this work. They must be posed some hypothesis, to get the properties a) b) and 
c) of the Poincare section and its return map F. 

The arguments above justify to wide the abstract mathematical model of a network of n inhibitory 
neurons, according to the following definition: 

Definition 2.15 The Abstract Mathematical Model. We say that a map F : B + B + , in a set 

B + homeomorphic to a compact ball of K™" 1 , models a generalized network of n inhibitory neurons if it 
verifies the statements a), b), c) of the Corollary |2.131 

3 The abstract dynamical system. 

Let B C M™ be a compact set, homeomorphic to a compact ball of R™ -1 . In particular Bi is connected. 
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Definition 3.1 A finite partition of B is a finite collection {-Bi}i<i< m of compact non empty sets Bi of B, 
homeomorphic to compact balls of M™ -1 , such that Ui<i<m Bi — B and int Bi n int Bj — 0, for i ^ j. 

Denote S = Ui^j Bi ^ j an d cau ^ the separation line, or Zme 0/ discontinuities , although it is not 
a line in the usual sense, but the union of the topological frontiers of Bi, each one homeomorphic to some 
(n — 2)-dimensional manifold. 

Definition 3.2 Given a finite partition {-Bi}i<j< m of B, we call F a piecewise continuous map on (B.V) 
with the separation property if F is a finite family F — {fi}i<i< m of homeomorphisms fi : Bi 1— > /i(-Bj) C B, 
such that fi(Bi) n fj(Bj) = if i ^ j. We note that F is multi-defined in the separation line S 1 . 
Each i?i shall be called a continuity piece of .F. 

Remark 3.3 A piecewise continuous map F with the separation property is globally injective because it 

is an homeomorphism in each continuity piece and two different continuities pieces have disjoint images. 

Therefore F^ 1 exists, uniquely defined in each point of F{B) = (J 4 fi(Bi). In fact: 

For any point x £ |J i fi(Bi), its backward first iterate is uniquely defined as F~ 1 (x) — (x), where i 

is the unique index value such that x S fi(Bi). 

Nevertheless F _1 is not necessarily injective because F is multidefined in S — ^ Bj)- 

F~ x is continuous in F(B), because F" 1 ^.^.) = /r and is an homeomorphism due to the Definition 

Definition 3.4 We say that F is uniformly locally contractive if there exists a constant < A < 1, called 
an uniform contraction rate for F, such that dist (fi(x), fi(y)) < A dist (x, y), for all x and y in the same 
Sj, for all 1 < i < m . 

Given a point x 6 _B, its image set is -F(a;) = {/i(a;) : ^ G -Bi}. If H C i?, its image set is F(H) = 
[j xeH F(x). We have that B D F(B) D . . . F k (B) D . . .. 

The second iterate of the point a; € B is the set i 7 ' 2 (a;) = F(F(x)). Analogously is defined the j— th. 
iterate as the set F J (x) for any j > 1. We convene to define = a; and F°(H) = H. 

Definition 3.5 For any natural number k > 1, we call atom of generation k to 

fik ° • • • ° /<2 ° 

where I = (ii,i2, ■ ■ ■ ,ik) S {1, 2, . . . , m} fe and Bi is the subset of Bi 1 where the composed function above is 
defined. (If B\ were an empty set, then the atom is empty.) Abusing of the notation we write the atom as: 

fi k °---°fi2° fii (Bh ) 

We note that each atom of generation k is a compact, not necessarily connected set, whose diameter is 
smaller than A^diamB. 

The set F k (B) is a compact set, formed by the union of all the not empty atoms of generation k. 
There are at most m k and at least m not empty atoms of generation k, where m is the number of 
continuity pieces of F. 

Definition 3.6 Given xq s B, a future orbit o + (xq) is a sequence of points {xi}i>o, starting in xq, such 
that Xi+i G F(xi) \fi > 0. Due to the multi-definition of F in the separation line S, the points of S and 
those that eventually fall in S may have more than one future orbit. 

A point y is in the limit set L + (o + (a:o)) of a future orbit of xq if there exists kj — > +00 such that 

xk 3 -> y- 

The limit set L + (xo) is the union of the limit sets of all its future orbits. 

The limit set L + (B) of the map F, also denoted as L + (F), is the union of the limit sets of all the points 

x e B. 
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Remark 3.7 Due to the compactness of the space B the limit set L + (o + (xq)) of any future orbit, is not 
empty. 

Also, it is standard to prove that L + (o + (xo)) is compact (because it is closed in the compact space B). 
Nevertheless L + [xq) may be not compact, if the point xo has infinitely many different future orbits. 
Finally, we assert that L + (o + (xq)) is invariant: L + (o + (xq)) ) = L + (o + (xo))- 

Proof: Consider y 6 L + (o + (xo)) . We have y = limj_>+oo Xk- € F{B) if kj > 1. 

F^ 1 : F(B) — > B is a continuous uniquely defined function in the compact set F(B) (see Remark |3 . 3|) . 
Then x kj -i = F~ 1 (x k] ) -» F~ 1 (y), so F^ 1 (y) e L+(o+(x )) proving that 

F- 1 (L+( +( a;o )))ci + ( + (xo)) 
Let us prove the converse inequality: F~ 1 ( L + (o + (x )) ) D L+(o + (xo)). 

i 7 = {fi : Bi i — ► _B} is defined and continuous in each of its finite number of pieces B t that are compact 
and cover B. Then there exists some i £ {1, 2, . . . , n} and a subsequence (that we still call Aj), such that 

y = lim Xfc, e Bi, Vj > : x kl e Bi, x kj +i = fi{x k .), fi(y) — lim/i(xfe,) = limxfc.+i 

There exists yi = /i(y) S -F(y) such that yi S i + (o + (xo)). In other words, j/ £ F^ 1 (L + (o + (xo)). This last 
assertion was proved for any y 6 L + (o + (a;o))- Therefore L + (o + (xq)) C -F _1 ( L + (o + (xo)) ) as wanted. □ 

Definition 3.8 We say that a point x is periodic of period p if there exists a first natural number p > 1 
such that x G F p (x). This is equivalent to x be a periodic point in the usual sense, for the uniquely defined 
map F~ l , i.e. F~ p (x) = x for some first natural number p > 1. 

We call the backward orbit of x (i.e. {^"^(x), j = 1, .. . ,p}), a periodic orbit with period p. 

We will prove in Lemma ?? that the limit set L + (B) is contained in the compact, totally disconnected 
set Kq = Hfc>i F k (B). It could be a Cantor set. But generically Kq shall be the union of a finite number 
of periodic orbits, as we shall prove in Theorem 14.11 

Definition 3.9 We say that F is finally periodic with period p if the limit set L + (F) is the union of only 
a finite number of periodic orbits with minimum common multiple of their periods equal to p. In this case 
we call limit cycles to the periodic orbits of F. 

We call basin of attraction of each limit cycle L to the set of points x 6 B whose limit set L + (x) is L. 



Topology in the space of piecewise continuous locally contractive maps in B. 

Let V — {-Bj}i<j< ro and Q = {^4i}i<i< TO be finite partitions (see Definition 13. ip of the compact region 
B with the same number m of pieces. 

We define the distance between V and Q as 

d{V,Q)= max Hdist(A,, Bi) (26) 

l<i<m 

where Hdist(A, B) denotes the Hausdorff distance between the two sets A and B. i.e. 

Hdist(A, B) = max{ dist (x,B), dist (y, A), x e A, y e B} 

Definition 3.10 Let F = {fi : Bi i— > 5}i<i< m and G = {gi : A.- L i— > S}i<,< m be locally contractive 
piecewise continuous maps on (B, V) and (B, Q) respectively. Given e > Owe say that G is a e-perturbation 
of F if 

max (gi- fdlBiHAi n o <e > |Af - A G | < e and d(P, Q) < e 

l<i<m 1 1 C 

where Af denotes t/ie uniform contraction rate of -F in its continuity pieces, defined in 13.41 and || • 1 1 c o 
denotes the C° distance in the functional space of continuous functions defined in a compact set K: 

IK.9 ~ J)\k He = max; dist (g(x),f{x)) 
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Definition 3.11 We say that the limit cycles of a finally periodic map F (see Definition 13. 9p are persistent 
if: 

For all e* > there exists e > such that all e-perturbations G of F are finally periodic with the same 
finite number of limit cycles (periodic orbits) than F, and such that each limit cycle Lq of G has the same 
period and is e*-near of some limit cycle Lf of F (i.e. the Hausdorff distance between Lg and Lf verifies 
Hdist(L G ,L F ) < e*). 

Definition 3.12 Denote S to the space of all the systems that are piecewise continuous with the separation 
property and locally contractive, according with the Definitions 13.21 and 13.41 

We say that a property P of the systems in S (for instance being finally periodic as will be shown in 
Theorem 14. lj) is (topologically) generic if it is verified, at least, by an open and dense subfamily of systems 
in the functional space S, with the topology (in S) defined in 13. 101 

Precisely, being generic means: 

1) The openness condition: For each piecewise continuous map F that verifies the property P there 
exist e > such that all e-perturbation of F also verifies P. 

2) The denseness condition: For each piecewise continuous map F that does not verify the property P, 
given e > 0, arbitrarily small, there exist some e-perturbation G of F such that G verifies the property P. 

The openness condition implies that the property P shall be robust under small perturbations of the 
system. It is robust under small changes, not only of a finite number of real parameters, but also of the 
functional parameter that defines the model itself. So the system should be structurally stable. When this 
robustness holds, the property P is still observed when the system, the model itself, does not stay exactly 
fixed, but is changed, even in some unknown fashion, remaining near the original one. 

The density condition combined with the openness condition, means that the only behavior that have 
chance to be observed under not exact experiments are those that verify the property P. In fact, if the 
system did not exhibit the property P, then some arbitrarily small change of it, would lead it to exhibit P 
robustly. 

The denseness condition implies that if the property P were generic, then the opposite property (Non-P) 
has null interior in the space of iS of systems, i.e. Non-P is not robust: some arbitrarily small change in 
the system will lead it to exhibit P. That is why we define the following: 

Definition 3.13 If the property P is generic, we say that any system that does not exhibit P is bifurcating, 
and Non-P is a not persistent property. 

4 The generic persistent periodic behavior. 

Theorem 4.1 Let F be a locally contractive piecewise continuous map with the separation property. Then 
generically F is finally periodic with persistent limit cycles. 

To prove Theorem 14. II we shall use the following lemma: 

Lemma 4.2 If there exists an integer k > 1 such that the compact set K — F k (B) does not intersect the 
separation line S of the partition into the continuity pieces of F, then F is finally periodic and its limit 
cycles are persistent. 

Proof: By hypothesis, dist (K, S) — d > 0, because K and S are disjoint compact sets. On the other 
hand K = F (B) = LUe.4 fc A where Ak denotes the family of all the atoms of generation k. 

As the diameter of each of the finite number of atoms of generation k is smaller than diam(B) X k , it 
converges to zero when k — » +oo. Thus, for all k large enough, it is smaller than d/2. 

We assert that each atom A of such generation k, is contained in the interior of some continuity piece 
Bi. In fact, fix a point x G A. As the continuities pieces cover the space B, there exists some (a priori 
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not necessarily unique) index i such that x G Bi. It is enough to prove that y G int(Bi) for all y G A 
(including x itself). 

We argue in the compact and connected metric space B, using known properties of any general compact 
and connected metric space, for instance the triangular property, and also the property asserting that the 
distance of a point y to a set, is the same that the distance of y to the frontier of that set. 

We denote Bf to the complement of Bi in B, and in the topology relative to B we denote: (Bf) to the 
closure of Bf, i.e the complement of int(Bi), and dBi to the frontier of Bi in B, dBi C S: 

dist (x, y) < diam(A) < d/2, dist (x,(Bf)) = dist (x,dB,)> dist (x,S)>d 

dist (y, (Bf)) > dist (x, (Bf)) - dist (x, y) > d - d/2 = d/2 > 

Therefore y g" (Bf) proving the assertion. 

We deduce that given an atom A G Ak, there exists and is unique a natural number such that 
A G int(Bi g ). Therefore F(A) is a single atom of generation k + 1. 

From the definition of atom in 13.51 we obtain that any atom of generation larger than k is contained 
in an atom of generation k. But each atom of generation k is in the interior of a piece of continuity of the 
partition {Bi}. 

We deduce that there exists a sequence of natural numbers {ih}h>0i such that 

A G mtB la , F(A) — f io (A) C intB n , F 2 (A) = f ix o f io (A) C intB is (27) 

and the successive images of the atom A of generation k, are single atoms of generation k + l,k + 2, . . . ,k + 
h, . . .. Therefore, the successive images of the atom A, in the sequence (f27|) . are contained in a sequence of 
atoms: A = Aq, A\, A2, . . ■ , Ah, ■ ■ ■ , all of generation k. 

The same property holds for any of these atoms of generation k, and each of them is contained in the 
interior of a continuity piece of F, so F is uniquely defined there and we have: 

A = A cintB io , F(A ) c Ax c intB H , F 2 (A Q ) C F(Ai) c A 2 C intB h , . . . , (28) 

The family of atoms of generation k is finite, so we conclude that there exists two first natural numbers 
< h< h +p such that F p (A h ) C A h . 

Note that, F p (Ah) is uniquely defined as fi h+p ° fi h+p _ 1 ■ ■ ■ fi h , because we are considering sets 
contained in the interior of the continuity pieces of F. 

Due to the uniform contractiveness of fi in each of its continuities piece, F p : Ah ► Ah, is uniformly 
contractive. The Brower Theorem of the Fixed Point states that in a complete metric space, any uniformly 
contractive map from a compact set to itself, has an unique fixed point, and all the orbits in the set converge 
to this fixed point in the future. Therefore, there exists in Ah a periodic point by F of period p > 1, and 
all the orbits with initial states in Ah have the periodic orbit L of p, as their limit set. 

By construction Ah was the image of A by an iterate F h uniquely defined. So we conclude that the 
limit set of all the points in the atom A is L. 

The construction above can be done starting with any initial atom A of generation k. And they are 
a finite family. We conclude that there exists one and at most a finite number of periodic limit cycles, 
attracting all the orbits of LUe.An, ^ = F k (B). 

The last assertion implies that the limit set of B is formed by that finite family of periodic limit cycles. 

Finally it is left to prove that the limit cycles are persistent according to the definition 13. Ill 

The condition of the hypothesis of this lemma is open in the topology defined in 13.101 because K and 
S are compact and at positive distance. 

We assert that the itinerary of each of the atoms A of generation k, for k fixed and large enough, remains 
unchanged when substituting F by G, being G a e-perturbation of F for e > small enough. 

In fact A and A (and also all the other atoms of generation k, with k fixed) are contained in the images 
by F k or by G k respectively, of some of their one-to-one corresponding continuity pieces. With k fixed, if 
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e > is sufficiently small, they remain at distance larger than the number 2d/3 > from the separation 
lines of F and of G respectively, being d = dist (K, S). 

Thus, if the generation k is chosen so the atoms have diameter smaller than d/3, repeating the argument 
at the beginning of this proof we show that A is in the interior of some continuity piece of F, and A is in 
the interior of the respective correspondent continuity piece of G. 

On the other hand, the future iterates of any atom of generation k by F, and also by G, are contained 
in the atoms of generation k. Therefore the images of an atom A or A of generation k, by all the future 
iterates of F or of G respectively, are in the interior of their respective one-to-one correspondent continuity 
pieces. Then the itineraries are the same as we asserted. 

As a consequence, the indexes io, i\, 12, ■ ■ ■ in the finite chain of atoms denoted in (|27| and (|28|) . remain 
unchanged, and therefore we deduce the following statement: 

A: The number of periodic orbits in the atoms of generation k, and their periods, remain unchanged, 
when substituting F by any e-perturbation G, if e > is sufficiently small. 

It is standard to prove by induction on k > 1 that for any e-perturbation G of F, such that A+e = A < 1, 
each atom A of generation k for G, is at distance smaller than Ylj^o^ 6 ^ < 2e/(l — A) = e* > of the 
respective atom A for F with the same itinerary. 

Therefore we deduce the following statement: 

B: Any periodic point found in an atom A of generation k for G, is at distance smaller than e* than 
the respective periodic point found in the correspondent atom A for F with the same itinerary. 

The statements A and B imply that the limit cycles are persistent according to Definition 13.111 □ 

Remark 4.3 In the proof of Lemma |4~2"1 we did not use the separation property fi(Bi)P\fj(Bj) = Vi ^ j. 
At the end of the proof of Lemma 14.21 we obtained that the piecewise continuous and locally contractive 
systems verifying the thesis of the Lemma I4.2[ even if they do not have the separation property, contain 
an open family of systems in the topology defined in 13.101 Then: 

In the space of all the piecewise continuous and locally contractive systems (even if they do not have the 
separation property), those whose limit set is formed by a finite number of persistent limit cycles form an 
open family. 

Nevertheless, to prove the genericity of the periodic persistent behavior, we need to prove that the 
family of periodic maps is dense in the space of systems. In the following proof, to obtain the density we 
shall restrict to the space of systems S that verify the separation property. 

Remark 4.4 From the proof of Lemma |4~2"1 the first integer k > 1 such that F k (B) f]S = may be very 
large, and so the period p may be very large. 

In fact, if the system has n « 10 12 neurons, and if no neuron becomes dead, i.e. it does not eventually 
remain forever under the threshold level without giving spikes, then the periodic sequences i±, . . . , i p , defined 
as the itinerary of the periodic limit cycles, have inside the period p, at least once each of all the indexes 
i G {1, 2, . . . , n}. Then p > n w 10 12 . 

As we have shown in the proof of the Lemma 12. 5[ there exists a minimum time T > between two 
consequent spikes. Suppose for instance that T ~ 10 [ms] and n « 10 12 . The lasting time of the periodic 
sequence could be approximately 10~ 3 x 10 12 [s] = 10 9 [s] > 31 years. So, if most of the neurons did not 
become dead, the observation of the theoretical periodic behavior of the inhibitory system in the future, 
could not be practical during a reasonable time of experimentation, and only the irregularities inside the 
period could be registered, showing the system as virtually chaotic. 

Proof of Theorem \4- 1[ Due to Lemma I4T21 the existence of a finite number of limit cycles attracting 
all the orbits of the space is verified at least for those systems in the hypothesis of 14.21 At the end of 
the proof of Lemma 14.21 we showed that its hypothesis is an open condition. To prove its genericity it is 
enough to prove now that the hypothesis of Lemma 14.21 is also a dense condition in the space of piecewise 
continuous contractive maps with the separation property. 
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Take F being not finally periodic. We shall prove that, for all e > there exists a e— perturbation G 
of F that verifies the hypothesis of Lemma [4~2l and thus G is finally periodic with persistent limit cycles. 
Let be given an arbitrarily small e > 0. 

The contractive homeomorphisms fi of the finite family F = {fi : Bi i— > B}i, with contraction rate 
< A < 1, can be C° extended to F e = {fi yt ■ Ui i— > B}i, where /i, e is an homeomorphism onto its image 
defined in compact neighborhoods Ui D Bi in B (i.e. Bi = Bi C int(Ui) C t/j = [/, C B where the closures 
and interiors of the sets are taken in the relative topology of B C K n ), such that fi^Bt — ft, and such that 
fi >e is still contractive, with a contraction rate < A' < 1 such that |A — A' < e. 

The extended map F e — {fi, e '■ Ui i— > -£>}i, is now multidefined on Ui/j ^ ^ ^ "S*- The separation 
property is an open condition, thus the extension F e still verifies /i, e (C/i) fl fj,e{Uj) = for all z ^ j, if the 
neighborhoods Ui and E/j are chosen at a sufficiently small Hausdorff distance from their respective pieces 
Bi and Bj, and e > is small enough. 

Call ei > to a positive real number smaller or equal than e, and also smaller or equal than the distance 
from Bi to the complement of Ui, for all i = 1, 2, . . . m. Precisely 



< ei = min{e, min dist (Bi, U?)} 

Ki<m 



Consider the compact sets: 



K + = f) (J fi k ,e0...of h>e (U h ) D K 

k>l (u,...,i fc )e{l,2...m} fc 

K=f] (J h k o...o f lx {B tl ) 

k>l (i u ...,i k )e{l,2...m} k 

Define the extended atoms of generation k > 1 for F e that form K + as fi ki€ o . . . o /t 1)£ (t/i 1 ) ) where 
(ii, 12, ■ ■ ■ , ik) is a word of length fc formed by symbols in {1, 2, . . . , m}. 

The diameter of each extended atom is smaller that diam(B) ■ -X' k . Therefore, for sufficiently large 
k > 1 all the extended atoms of generation k that form K + have diameters smaller that ei/2. 

We assert that the extended atoms of generation k > 1 are pairwise disjoint: in fact, for two different 
i ^ j the images are disjoint: fi, e (Ui) fl fj,e(Uj) = 0. So the atoms of generation 1 are pairwise disjoint. 
Two extended atoms of generation k are fi k<e o . . . o f ix £ (JJ ix ) and /j fc !e o . . . o fj ± _ e (Uj 1 ) . They can intersect 
if and only if (i\,i2, ■ ■ ■ ,ik) — {ji,j2, ■ ■ ■ ,jk) because each fi t£ is an homeomorphism onto its image. So, 
they intersect if and only if they coincide. 

By construction, Ui D Bi and fi, e \Bi = fi- Therefore each of the atoms of generation k for F, is contained 
in the respective extended atom of generation fe for F e , that has the same finite word (i\, 12, ■ ■ ■ , ik)- 

If none of the extended atoms of generation fc intersects S, then none of the atoms of generation fc for 
F intersects S, and the system verifies the hypothesis of Lemma 14.21 So, in this case, there is nothing to 
prove, because the given system F verifies the thesis of the Lemma B~2l and thus, it is finally periodic with 
persistent limit cycles. 

If some of the extended atoms of generation k intersects S, consider a new finite partition Q = {Ci}i<i< TO 
of B such that the distance, defined in ([2^]) . between Q and the given partition V of F, is smaller than 
ei > 0: 

dist (V, Q) < ei < e 

We shall besides choose the new partition Q such that the new separation line Sq = Ui^ (Ci ^ Q) 
does not intersect the extended atoms of generation k of K + . This last condition is possible because the 
diameters of the generalized atoms are all smaller than ei/2, they are compact pairwise disjoint sets, and 
the distance between the two partitions V and Q (which is smaller than t\ > 0) can be chosen larger than 
ei/2, defined in (|26p as the maximum Hausdorff distance between their respective pieces. (We note that 
the old, and principally the new, separation lines S-p and Sq, are not necessarily C 1 nor even Lipschitz 
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manifolds in the space B, and even if they are, they do not need to be C 1 or Lipschitz near one from the 
other, to be near with the Hausdorff distance). 

The first condition dist (V, Q) < e\, joined with the assumption dist (Uf,Bi) > e\, where Bi is the 
i — th piece of the partition V , implies that the respective piece Ci of the partition Q verifies Ci C Ui. 
Therefore the extension fi. t : Ui i— > B in F e can be restricted to Ci. 

Define G — {gi : Ci i— > £?}i< m where <?.; — /i, e |cv By construction G and F coincide in Ci D Bi, the 
distance between the respective partitions V and Q is smaller than ei < e, and the difference of their 
respective contraction rates A' and A is also smaller than e. So G is a e-perturbation of the given F, 
according to the Definition 13.101 It is enough to prove now that G is finally periodic with persistent limit 
cycles. 

Consider the limit set Kq of G as follows: 



As G is a restriction of F e to the sets Ci C Ui, we have that Kg C K + , and in particular for all k > 1 
the atoms of generation k for G, i.e. gi k o . . . o (C,^), are contained in the extended atoms of generation 
k for F e . 

By construction the separation line Sg = Sq among the continuity pieces Ci of G is disjoint with the 
extended atoms of generation k of F t . Therefore, it is also disjoint with the atoms of generation k of G. 
Then G k (B) f] Sg — and, applying lemma ET~2l G is finally periodic with persistent limit cycles. □ 

5 Open mathematical questions. 

It is possible (but not immediate) to construct, in a compact ball B of any dimension n — 1 > 2, piecewise 
continuous systems, uniformly locally contractive and with the separation property, as defined in Section [3j 
that do not verify the thesis of the Theorem 14.11 and thus their limit set is not composed only by periodic 
limit cycles. 

Suppose that the system had C r regularity, with r = Lipschitz or with r > 1, i.e. the continuity pieces 
Bi and the separation lines that form S, are bi-Lipschitz homcomorphic (or C r diffcomorphic respectively) 
to (n — 1) or (n — 2) dimensional balls or manifolds, and the homeomorphisms fi in their continuity pieces 
Bi, are bi-Lipschitz (or C"*-diffeomorphisms respectively). 

With this additional assumption of regularity of F, it is an open question to construct examples that 
do not verify the thesis of Theorem 14. 11 In other words, assuming more regularity, it is unknown if the 
system has to exhibit a limit set always formed only by periodic orbits. 

On the other hand, it is also unknown if the existence of a finite number of limit cycles, as in Theorem 
14.11 is generic for Lipschitz or C 1 regular, locally contractive and piecewise continuous systems with the 
separation property. 

6 Conclusions 

The discontinuities of the Poincare transformation F, due to spike phenomena in the neural network, play 
an essential role to study these systems, although it is an obstruction to apply mostly previously known 
results of the Theory of Dynamics Systems, which is mostly developed for continuous dynamics. 

In the generic stable case, the recurrent orbits are all periodic, and all the initial states lead to limit 
cycles. 

Due to the non-genericity of the bifurcating case, which is a consequence of Theorem 14 .11 those dynamic 
would never be seen in experiments: in fact, arbitrarily small perturbations in the parameters of the system 
will lead it to periodic or quasi- periodic dynamics. These perturbations stabilize the system, to exhibit a 
limit set composed only by limit cycles. 






k>l (i 1 ,...,i fe )e{l,2...m}* 
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Also we showed that the inter-spike interval is bounded away from zero. It means that, generically, 
when the system is periodic, in spite of having preferred periodic patrons of discharges, the neurons do not 
synchronize in phase. 

On the other hand, if the number of neurons in the system is very large, the limit cycles of the network 
may have a very large period p, much larger than the observation time, or even than the life time of the 
biological system. Therefore, in spite of being asymptotically periodic, these systems may never show its 
regularity. These two facts: extremely large periods, and irregularity inside the period, allow us to assert 
that those persistent systems with very large period p shall be in fact non-predictible for the experimenter, 
and will be perceived as virtually chaotic. 
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